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The existence of a primordial magnetic field (PMF) would affect both the temperature and po- 
larization anisotropics of the cosmic microwave background (CMB). It also provides a plausible 
explanation for the possible disparity between observations and theoretical fits to the CMB power 
' spectrum. Here we report on calculations of not only the numerical CMB power spectrum from the 

I PMF, but also the correlations between the CMB power spectrum from the PMF and the primary 

curvature perturbations. We then deduce a precise estimate of the PMF effect on all modes of 
' perturbations. We find that the PMF affects not only the CMB TT and TE modes on small angular 

^ ' scales, but also on large angular scales. The introduction of a PMF leads to a better fit to the CMB 

power spectrum for the higher multipoles, and the fit at lowest multipoles can be used to constrain 
\^ • the correlation of the PMF with the density fluctuations for large negative values of the spectral 

index. Our prediction for the BB mode for a PMF average field strength \B\\ = 4.0 nC is consistent 
with the upper limit on the BB mode deduced from the latest CMB observations. We find that the 
BB mode is dominated by the vector mode of the PMF for higher multipoles. We also show that 
by fitting the complete power spectrum one can break the degeneracy between the PMF amplitude 
and its power spectral index. 
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■ PACS numbers: 98.62.En,98.70.Vc 

^ ■ I. INTRODUCTION 

> 

\ Magnetic fields in clusters of galaxies have been observed P, d, H, Q with a strength of 0.1 — 1.0 ji G. The existence 
lO • of a primordial magnetic field (PMF) of order 1 nG whose field lines collapse as structure forms is one possible 
\ explanation for such magnetic fields in galactic clusters. The origin and detection of the PMF is, hence, a subject of 
considerable interest in modern cosmology. Moreover, the PMF could influence a variety of p henomena in the early 
universe [| such as the cosmic microwave background (CMB) [1, 0, i, i, [13, El, [H, [II, III [H, [H, 113, El , or the 
OO ; matter density field [H, HO, IMI, HI . 

' Temperature and polarization anisotropies in the CMB provide very precise information on t he p hysical processes 
' ' ■ in operation during the early universe (WMAP [2l,|2l,|2|, ACBAR [23, CBI [Hill, DASI HO], BOOMERANG 
[Slj , and VSA [33| ) . The CMB power spectrum from ACBAR and CBI, has indicated a potential discrepancy between 
. , these observations at higher multipoles ^ > 2000 and the best-fit cosmological model to the WMAP power spectrum. 

■ A straightforward extension of the fit [2^ to the WMAP data predicts a rapidly declining power spectrum in the 
■ ■ ' large multipole range due to the finite thickness of the photon last scattering surface and the Silk damping effect. 

The ACBAR and CBI experiments, however, indicate continued power up to ^ ~ 4000. This dis crep ancy is difficult 
to account for by a simple retuning of cosmological parameters or by the Sunyev-Zeldovich effect |23l. Issl. |33| . Among 
other possible explanations, an inhomogeneous cosmological magnetic field generated before the CMB last-scattering 
epoch provides a plausible mechanism [35j to produce excess power at high multipoles. Such a field excites an Alfven- 
wave mode in the primordial baryon-photon plasma and induces small rotational velocity perturbations. Since this 
mode can survive on scales below those at which Silk damping occurs during recombination [13l . |36| , it could be a new 
source of the CMB anisotropies on small angular scales. The present work, therefore, is an attempt to more precisely 
study the evolution of cosmological perturbations with a PMF. 

Previous work [1, 0, i, i, E, El, 111, El, El, El, El, El, El, IH, IH, HI has shown that one can obtain information 
about the PMF from the CMB temperature anisotropies and polarization. However, in those works attention was only 
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given to a subset of the modes of the CMB anisotropies. In the present work, therefore, we study the comprehensive 
effect of the PMF on all modes of the CMB perturbations. Furthermore, in order to clarify the role of the PMF in 
the CMB, we take into consideration the possible correlation between the CMB fluctuations induced by the PMF and 
those due to primordial curvature and tensor perturbations. Also, we numerically evaluate the CMB power spectrum 
from the stochastic PMF and thereby avoid recourse to analytic approximations. 

In this article, we use adiabatic initial conditions for the evolution of primary density perturbations and consider 
isocurvature (isothermal) initial conditions when estimating effects on the CMB anisotropy induced by the PMF 
11 , . Throughout this article we fix the best fit cosmological parameters of the ACDM + Tensor model as follows 
H: h = 0.792, ^IbK^ = 0.02336, ^mh? = 0.1189, ns = 0.987, r = 0.55, ut = -r/8 = -0.069, and = 0.091 in flat 
universe models, where h denotes the Hubble parameter in units of 100 km s^-'^Mpc""'^, f2f, and fJ™ are the baryon and 
cold dark matter densities in units of the critical density, ng is the spectral index of the primordial scalar fluctuations, 
r the ratio of the amplitude of the tensor fluctuations to the scalar potential fluctuations, tit is the spectral index of 
the primordial tensor fluctuations, and Tc is the optical depth for Compton scattering. 



II. PRIMORDIAL MAGNETIC FIELD 



Before recombination, Thomson scattering between photons and electrons along with Coulomb interactions between 
electrons and baryons were sufRciently rapid that the photon-baryon system behaved as a single tightly coupled fluid. 
Since the trajectory of plasma particles is bent by Lorentz forces in a magnetic field, photons are indirectly influenced 
by the magnetic field through Thomson scattering. Let us consider the PMF created at some moment during the 
radiation-dominated epoch. The energy density of the magnetic field can be treated as a first order perturbation 
upon a flat Friedmann-Robertson- Walker (FRW) background metric. In the linear approximation, the magnetic field 
evolves as a stiff source. Therefore, we can discard all back reactions from the magnetohydrodynamic (MHD) fluid 
onto the field itself. 



A. Power Spectrum from the PMF 

To derive the power spectrum from the PMF we begin with the electromagnetic tensor in the usual form 
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where Ei and Bi are the electric and magnetic fields. [Here we use natural units, i.e. c = h = 1.] The energy 
momentum tensor for electromagnetism is 

^"^[EM] - ^ (^F'^-'F^^ - Ig'^I^F^sF'''^ . (2) 

The Maxwell stress tensor, a^^, is derived from the space-space components of the electromagnetic energy momentum 
tensor, 

Tik ik 
[EM] = cr — 

4— Ie'e'' + B'B'' - -y^iE^ + B'^)\. (3) 

a'' Att I 2 J 

As mentioned above, within the linear approximation [s^ we can discard the MHD back reaction onto the field 
itself ^ij. The conductivity of the primordial plasma is very large, and is "frozen- in" [l^, This is a very good 
approximation during the epochs of interest here. Furthermore, we can neglect the electric field, i.e. ~ 0, and can 
decouple the time evolution of the magnetic field from its spatial dependence, i.e. B(t, x) — B(x)/a^ for very large 
scales fllj. In this way we obtain the following equations, 

B^ 

^°°[™1 = 8^ ' ^"^^ 
T'°[EM] = T'-"'[em] = , (5) 
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-T'\EM]^a''^-^{2B'B''-y''B^) . (6) 
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We assume that the PMF Bq is statistically homogeneous, isotropic and random. For such a magnetic field, the 
fluctuation power spectrum can be taken as a power-law S{k) =< B{k)B* (k) >(x A:"^ [T^ where rie is the power-law 
spectral index of the PMF. The index n-Q can be either negative or positive depending upon the physical processes of 
magnetic field creation. From Ref. [l^l , a two-point correlation function for the PMF can be defined by 

B\k)B' ik')) = ^-A___^k--p^^^k)6ik-k% k<kc , (7) 



where 



P'\k) = S'^--^ . (8) 



Here, B\ is the magnetic comoving mean-field amplitude obtained by smoothing over a Gaussian sphere of comoving 
radius A, and k\ = 27r/A (A = 1 Mpc in this paper). The cutoff wave number kc in the magnetic power spectrum is 
defined by 



kfj ^("^dcc), ^ Tdcc, 

where is the mean free path of photons, and rdoc is the conformal time of the decoupling of photons from baryons. 

B. Scalar Mode 

We obtain the power spectrum of the PMF energy density |i?[EM:S](^j ^ ^) a-nd the Lorenz force for the 
scalar modes \I{^^yi.'g,-^{k,T)\'^5{k— k') according to the following relations, 

|%M:S](fc,T)|25(fc-fc')= (^(7^(fc,^)[EM:Sl]r*(fc',T)[EM:Sl]> , (10) 

and 

|n[EM:S](fc,T)|^<^(fc- k') = J^^{{Tik,T)[-EM:Sl] " T {k, t) [EM:S2]) 

X (T*(fc',T)[EM:Sl] -r*(fc',r)[EM:S2])) • 

(11) 

In the case of a power law stochastic magnetic field, an explicit expression to evaluate the ensemble averages for the 
above spectra is given by 

mi. ^ T*( ^ \ 1 [ (2^)"B+8 Bl Y 

{1 [k, r)[EM:Si] J (P, T)[EM:Si]) - 24(27r)8a8 1 '^k^+^rJn^M) J 



kr^-'''' (t) = I i^p+p) lo dr'^, r < tjcc (g) 



.JdW^^\k^k'r[i + i^;J0f]sik^p) 



(12) 



Equivalently, in terms of an angular integration one can write, 

if (27r)"=^+® B"^ Y 
{T{k, t)[em:Si]T{p, r)fEM:Si]) = 2^2TTya» \ ~2k^^f~(Ji^ J 

X J dk'k"'^+^ J dClk-k'l"^-"^ {{1+C^)e ~4kk'C + 2k'^}S{k-p) 



(13) 



where we define 



i i' k'k 

C = cose = fe • = — — . (14) 
k'k 
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In almost all previous work the sum of the terms in brackets which include C in the k integral of Eq. ((T3l) have been 
set to unity. In this paper, however, we calculate Eq. (|13p explicitly by integrating all of the terms. In this way we 
obtain the following expression. 



liB+5 



1 f (27r)": 

X / ^fc''^'"^^' i^'Z^' {{k + k') 

J [fcfc'nB(nB + 2)(nB +4) ' 
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(15) 



A similar calculation gives the power spectrum of the PMF tension and the power spectrum of the correlation between 
pressure and tension as follows, 

(%M:S2](fc)T[*EM:S2](^)) = ^kl''"''^ 



and 
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,,,,,, , ^ f {(fc + ^■')""+' -\k- A:'|"-+n 

(fcA:')^'^B(?^B + 2)(nB +4) L ' ' ' ^ 

(7'[EM:Sl](fc)7'[*EM:S2](^)) + (%M:S2](fc)r[EM:Sl](^)) 

1 /(27r)' 



V'B+5 ^2 ^ ^ 



X / dk'k"'^+^ 



^|(fc + fc')"B+3_ |^_^,|nB+3} 

_(fc/^')^"B(?^B + 2) ^ ' 'J 

"Pfc'3nB(nB+2)(nB+4) ^ ~ + + ^'^"""'^ 

'/fc3fc'2„B(nB+2)(nB+4) + " " ^''"^'^ 

^/fc3fc'4nB(nB+2)(nB+4)(nB + 6)^^^^ ' ' > 



(16) 
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C. Vector Mode 

We obtain the power spectrum of the PMF Lorenz force for the vector mode from the following 

1 



|n[EM:V](fc)P<5(fc-fe') 



2(27r)- 



{Ti[EM:V]{k)T* [EM:V]{k')) 



(18) 
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For the case of a power law stochastic magnetic field, an explicit expression for the ensemble average is given by 

[EM:V] 

1 /(27r)"=^+8 Bl 



4(27r)7a8 2k 



"3+3 p / Kb +3 'I 
A ^ V 2 ; 



J dk' J dCk'^'^+'^ik- k')"""-"^ 



(1 - C^) {2fc2 - bCkk' + {2C^ + l)fc'2} . 



Integrating Eq. (fT9|) over C, after a lengthy calculation, we obtain the following expression. 
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+ (llfc'2 + l^kk' + 2fc2)(/j + /j')«B+4} 
1 
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^(24fc'2 + 15A:/c')(A: + /c')"''"^^} 
24 
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D. Tensor Mode 



We obtain the power spectrum of the PMF for the tensor mode according to the following 



|n[EM:T]|''5(fc-fc') 
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4(27r); 



(r'^(fc)[EM:T]7ji*(fc')[EM:T]) ■ 



(21) 



For the case of a power law stochastic magnetic field, we use the explicit expression for the ensemble average as given 



(T^^ ik)[-EM:T]Tij* (k) [EM:T] ) 



(27r)"B+8 ^2 



4(27r)7a8 [ 2^5^^+^ ^ i'-^) 



dk'k"'''\k- k'\''^{l+C'} { 1 



/\2 



{k^Ck') 
Ifc- fc'P 



(22) 



6 



Again integrating Eq. (|22p over C. after a lengthy calculation, we obtain the following, 
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(23) 



These relations (Eqs. [TU] - 1^5)) constitute the various components of the PMF power spectrum evaluated in the 
present work. 



E. Numerical vs. Analytical 

In previous analyses @, [13, [O, [13, 07 SI] an analytic approximation to the power spectrum of the PMF was utilized. 
Specifically, {k + k'Y^+'^ -\k- k'\'^^+'^ ~ 2(nB + 2)fc"=^+iA:' for k' < k, and 2(nB + 2)fcA:'"=^+i for k' > k. However, this 
is only a good approximation for fc' ^ fc or fc' 3> k. Outside of this range it can be a poor approximation since the 
neglected terms are not small in general. In previous work the source power spectrum from the PMF was obtained 
by integrating the approximate equations over < k' < kc, thus leading to smaller values than those evaluated 
numerically for tib < -2 as displayed in Figs.[T^ and[T|D. For example, when ub = —2.9, (and k < 0.1) the relative 
errors in the squared power spectrum from the PMF (iiimM/^APp) can be as large as 400%, 20%, and 500% for 
the scalar, vector and tensor modes, respectively. Figure lb, shows the same result for the more familiar multipole 
coefficients, Ci, which are related to the various components of the power spectrum via an integration over wave 
number, 

{21 + lfcf'\PMF) = 1 1 dkk'[@^^^^{k)e^^^Uk)] • (24) 
Here, O is the photon moment, 

e = T(fc) X (H(fc)2)V2 ^ (25) 

where T{k) is the transfer function, and the index X = S, V, or T denotes the scalar, vector, or tensor modes, 
respectively. Lines are drawn on Figure lb for ^ = 10 and 1000, which roughly corresponds to the k = 0.1 Mpc^^ and 
0.001 Mpc~^ shown on Figure la. 

Figures la and lb show that the deviations of the previous approximation from our numerical estimation increases 
as riB decreases from —1.5 to —3.0. The deviations of the tensor and scalar modes in particular arc as much as one 
order of magnitude near ub = —3.0. We understand this for the following reasons: Because integrands of all modes 
are dominated by values for range of fc' < fc for ub < —1.5, we can consider only this range. Since k and k' are 
smaller than unity, and for ub < —2, values of \k — k'\"^~^'^ around k ^ k' dominate the power spectrum of all modes 
fEas. ll5lTBlll7[E5)) except the vector mode (Eq. EO)) . the deviations of the previous approximation from our numerical 
estimation increase exponentially for tt-b < — 2 (Fig. [T^b) [63|- Since the order of the PMF source power spectra of 
all modes ([121 [H [13 [201 USD is 2nB + 3, the values for k' > k and riB > -1.5 in integrands of all modes (fT51 [TBI [T71 
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[20l [23]) dominate the result of these integrations. Therefore, if the integral range extends beyond k for ub > —1.5, we 
can ignore the integrands of all modes for small enough k' and the previous approximation is adequate for the PMF 
source power spectrum. Thus, ratios of all modes are constant for tt-b > —1.5. These results are almost the same as 
in Brown and Crittenden [46t [6l|. The constant ratios of all modes for hb > —1.5, however, are not unity because 
the terms including cosine factors [e.g. sum of terms within the bracket of Eq. (|15p ]. is not unity as was assumed in 
the previous approximation. 

There is, however, only a slight deviation for the vector mode (cf. the middle panels in Figs. [1^ and[T|3). This 
is because the power-law index of the k — k' term in the PMF source power spectrum of the vector mode scales as 
(cf. Eq. [20|l riB + 3; while the power-law index for the {k — k') term for the scalar and tensor modes is ne + 2 in Egs.fTSl 
and [231 Also, the k + k' term slightly dominates the PMF source power spectrum of the vector mode (n[EM:V]{l^)) 
for hb < —1.5. Hence, unlike the scalar and tensor modes, there is only a slight deviation in the vector mode for 
riB < —2. Additionally, since k and k' are smaller than unity and the power-law index of the (fc -I- k') term, which 
dominates n[EM:V](^ = 0.1), is positive, IIjemiV]!^ = 0-1) > n[EM:V](^ = 0.001). This is in contrast to the scalar and 
tensor modes for which lliEM:V]{k = 0.1) < lliEM:V]{k = 0.001). The reason for this can be traced to the negative 
value of the power-law index (ub + 2) for the scalar and tensor terms when hb < —2. 

Also note on Figs.[T^ and[T|3 that even in the limit of riB > —2, the deviations to not asymptotically approach unity. 
The reason for this is due to the inclusion of the correct C factors in the numerical angular integrals for Eqs. (fT3l [T9l 
and l22p which are not unity as is assumed in the analytic approximations. 



III. EVOLUTION EQUATIONS 



In this work, we have implemented a numerical method ([22, |4^) to evaluate the PMF source power spectrum. 
Using this method, we are able to quantitatively evolve the cut off scale and thereby reliably calculate the effects of 
the PMF on the observed CMB power spectrum. We now summarize the essential evolution equations for each mode. 



A. Scalar Mode 

For the scalar mode we obtain the following equations in fc-space [l^, lU I3 IS US [HI : 



k^(j) + 3H{^ + Hi;) = A7:Ga^ {E[EM:S]{k,T) - dptot} (26) 
k\<j)-^) = -127rGa2|Z[EM:S](fc,T)-(p, + P,)a,} 

= -UnGa^ I^^E[EM:S]{k,T) +n[EM:S]{k,T) - {p, + P,)<j,^ (27) 
<5(S) = _(i + ^) (y(s) + - 3^ - ^) -5^'^ 

-^{E[EM:S]{k,T)+6HE[EM:S]{k,T)^ ' ^^^^ 

= _i/(i_3^)„(S)„^„(S) + ^^!^_fc2^^^2^ 
1 + w op I + w 

+fc2HM^^ , (29) 
47r/9 

where w = p/p. Note that for the photon (S^^' = 48g^'', and v^'' = fcO^^'. Massless neutrinos obey Eqs. ([28)1 and 
([^^l without the Thomson coupling term. In the continuity and Euler relations (Eqs. [551 and [23) for the scalar mode, 
we can just add the energy density and pressure of the PMF to the energy density and pressure of cosmic fluids, 
respectively. Since the baryon fluid behaves like a nonrelativistic fluid during the epoch of interest, we may neglect 
w and SP^^^ /6p'j^\ except the acoustic term Csk'^5f'\ Also, the shear stress of baryons is negligible Since we 

concentrate on scalar type perturbations in this paper, we do not consider the magneto-rotational instability from 
the shear stress of the PMF and baryon fluid [l^l . 

From equations ([28)) and ([29]). we obtain the same form for the evolution equations of photons and baryons as in 
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previous work [13, IH, |4^, [s^, [5l| , by considering the Compton interaction between baryons and photons. 



e(b + 3H{<P + Hi,) = AnGa^ {E[EM:S]ik,T)^Sp,,,} (30) 

= -UirGa^ I^^E[EM:S]{h t) +n[EM:S]{k, t) - {p, + P,)a,^ (31) 

Cm = -«g^M+30, (32) 

<5f = -^t'f +40, (34) 

-5^^) = -^«F+40, (35) 

- fc^ Q<5(S) _ a,) + an,aT(^;f - «f )) + fc^V, (36) 

= fc^Q5(S)-a.)+fc2V, (37) 

^ _^;(S)+30 (38) 

■(S) « (S) , 2;2r(S) , 4p^ (gj (s), 2; 

a 6ph 

, 3 2n[EM^S](fc^ , . 



where R = (3/4)(ph/p^) is the inertial density ratio between baryons and photons, Ue is the free electron density, ar 
is the Thomson scattering cross section, and cr^ of the second term on the right hand side of equation ([36|) is the shear 
stress of the photons with the PMF. Since ne < is favored by constraints from the gravitational wave background 
p^ l and the effect of the PMF is not influenced by the time evolution of the cut off scale kc for this range of ub , we 
approximately set i?[EM:S] a~'^ in the following analysis. 



B. Vector Mode 



The evolution of the vector potential F(t, k) under the influence of a stochastic PMF can be written [43, [4^ as 

leTra^GlIfEMiV] (k, t) 



V + 2-V 
a 



k 

^^^^2Pj]h+P^ (40) 



where the dot denotes a conformal time derivative, while pi and tt, are the pressure and the anisotropic stress of the 
photons (i = 7) and neutrinos (i = v). Here, we have omitted the vector anisotropic stress of the plasma which is 
negligible in general. In the absence of a magnetic source term, the homogeneous solution of Eq. (|40|) behaves like 
V (X We take a oc t during the radiation-dominated epoch. The magnetic field, therefore, causes the vector 

perturbations to decay less rapidly (oc 1/a instead of 1/a^) with the universal expansion. 

Since the vector perturbations cannot generate density perturbations, we have = Sj^^ ~ 0, where S'^^ and sj^^ 
are the perturbations of the photon and baryon energy densities, respectively. 

The magnetic field affects the photon-baryon fluid dynamics via a Lorentz force term in the baryon Euler equations. 
Following [13, SI], the Euler equations for the neutrino, photon and baryon velocities, , and v^'^ are written 
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as 



e'-v^^-M^e'T) , (41) 



a , m 1 



T/ + f.(«(v)-.rO--M^©V2M > (42) 



(V) -,> , "/ (V) -L . . (V) (V)n 

b - ^ + -(«b - ^) - ;^^=(^7 ^ «b ) 

3 n[EM:V](k, t) 
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(43) 



For the photons w^^^ = 0i^\ while O^J^^ and 6^^^^ are quadrupole moments of the neutrino and photon angular 



distributions, respectively. These quantities are proportional to the anisotropic stress tensors. Equations 141 
denote the vector equations of motion for the cosmic fluid, which arise from the conservation of energy-momentum. 

C. Tensor mode 

The Einstein equations tell us that the tensor mode H is governed by [13, |H, |4^, [5l| 

H + 2^H + k^H - SnGa^ (jIiem-.t] + ^''P) , (44) 

(T) 

where III is the anisotropic stress for neutrinos. 

D. Initial Conditions 

We need to specify the initial perturbations for solving the evolution equations derived in the previous sections. 
We start solving the evolution equations at early times when a given k mode is still outside the horizon, i.e. the 
dimensionless parameter fcr <JC 1. We consider only the radiation-dominated epoch since the numerical integration 
for all of the k modes of interest will start in this era. Baryons and photons are tightly coupled at this early time. 
The expansion rate is H = . We then derive initial conditions for all of the modes utilizing the method of 
Refs. [ll|,[2i,[43,[4i,[4i[53,[5l|. 

1. scalar mode initial condition 



We can assume that all density fields are zero since initially, since the PMF only affects the velocity field of ionized 
baryons, i.e. the Lorenz force, and the density fields are not directly affected by the PMF. In the radiation dominated 
era, photons and neutrinos are important in the RHS of Einstein equations. Equations for photons and neutrinos are 

= -^«f -f4(^, (45) 

S^J"^ = -lv(J^^+4, (46) 



(47) 



€^ - (jS^J^^ - a^f^ + fcV) , (48) 
cri'' = ^-f^. (49) 
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Here we omitted higher multipolc moments £ > 1 for photons and i > 2 for neutrinos. At the lowest order in kr, 
initial conditions of Eqs. (|30ll551) 



where 



CTB 

Rb 



^[EM:S] 

Pi 



Pj + Pv 

Pv 
P-i + Pu' 
-£'[EM:S] 

Pi 



2. vector mode initial condition 
At early times we neglect the vector anisotropic stress of the plasma, which is in general small. Equation (|40p gives 

y(r,k) + 2V(r,k)^-^'"°'^"f"-^'^"'^^ , (55) 
a k 

which can be easily solved to obtain [l^ 



167rG'a^n[EM:V](T,k) _ ^^^^ 



From Eqs.dUl) and (gS]), we obtain [ll|,ll 



(V) _ 3 n[EM:V](T, k)fcT ^ leTrGa^njEMiV] (t, k) ^ 
~ 4 p^(l + i?) A: 



(57) 



3. tensor mode initial condition 
From Ref . [lli , we obtain the initial condition of the tensor mode as follows 



^ 3R-, In ^— ) (58) 

V'''PMF/ 

15 R^ \ n[EM:T] .2 



n = 1 + - tl}E±Ii^krf (59) 

^ V 8i?^ + 30y 28 4i?^ + 15 J ^ ^ ^ ' 

^^i' - 1 o .p , ,K - 1 - , -1^ -5- f — n[EM:T] , (60) 



3 4i?^ + 15 V 14 4i?^ + 15/ R^j f™^^ 
where r"^ is the time of neutrino decoupling, and tpmf is the time of generating the PMF, R^ = fl-f/^Q^, + fi^), and 

R^ = n^/{fi^ + n^). 
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IV. CORRELATIONS IN THE POWER SPECTRA 

Although possible origins of the PMF have been studied by many authors, there is no consensus as to how the PMF 
correlates with the primordial density fluctuations. Nonetheless, almost all previous works have assumed that there is 
no correlation between them In order to study the possible effects of a PMF in a more general manner, we here 
consider possible correlations between the PMF and the primordial density and tensor fluctuations. To do this we 
introduce a coefficient s^^-* to parameterize the correlation between the PMF source and the primary power spectrum 
[13, E^l- The generalized multipole coefflcients c/'^'' (with X = S or T for scalar or tensor modes) then become, 

+.(^){ef )(fc)e[^S*p(fc) + e[^Sp(fc)ef (ei) 

where Q is the photon moment as deflned in Eq. 1251 The last term is the correlation between the primary fluctuations 
and the power spectrum from the PMF. Among the many possible sets of two correlation coefficients that satisfy 
— 1 < s^^"^ < 1, the two cases of s^^^ = s^"""^ = ±1 are expected to show the maximum absolute effects from the PMF. 
These limits represent the effective range of the resultant CMB anisotropics when a PMF is present. 

V. RESULTS AND DISCUSSION 

We have explored effects of a PMF on the CMB for the allowed PMF parameters which were deduced in our 
previous work [l3| (i.e. B\ < 10 nG and ?ib < —2.4) . The upper panel of Figure 2 illustrates the CMB temperature 
and polarization anisotropics from the PMF for scalar, vector and tensor modes for the case when B\ = 4.0 nG and 
riR = —2.9 or —2.5 as labeled. The scalar mode dominates for lower £ of the TT and TE modes as shown by Giovannini 
[40]. In particular, it is comparable in power to the primary TT mode for {Bx,nB) = (4.0 nG,— 2.9). We note that 
the curves with (Ba, «b, s^^\ s^"""-*) = (4.0 nG, —2.5, 1, 1) give the best fits to the observed power spectra in both the 
re gion s of high and low £. This result is complementary to and consistent with the nucleosynthesis constraints derived 
in (i^l as shown in [4^ . 

For illustration, let us assume the same scale-invariant power spectrum for both the PMF and the primordial 
curvature perturbations. In this case, the ratio of the density and velocity perturbations induced by the primordial 
curvature perturbation to those by the PMF is proportional to k^. Therefore, the temperature anisotropics from the 
PMF are larger for lower £ compared with those from primordial curvature perturbations. Furthermore, the power 
of the CMB temperature anisotropics from the PMF for lower £ depends not only on B\ but also strongly on ??b 

[Panel (la) of Fig. 2]. Note, that the magnetic field which is continually sourcing fluctuations docs not spoil the phase 

(s) 

coherence (cf. cosmic defects). The basic behavior of acoustic oscillations is affected by the pressure of the fluid kcsSl 
and the potential kcj). Since the pressure of the PMF is sufficiently less than the thermal fluid pressure, the PMF dose 
not affect the phase coherence signiflcantly. 

The PMF also affects the CMB power spectrum on small angular scales for two reasons. First, the PMF energy 
density fluctuations depend only on the scale factor a and can survive below Silk damping scale. Therefore, the 
PMF continues to source the fluctuations through the Lorcntz force even below the Silk damping scale. Second, the 
vector mode from the PMF can be larger than the scalar modes both from the PMF and primordial perturbations 
at small scales. This is because, after horizon crossing, the latter cannot grow due to the photon pressure leading 
to acoustic oscillations, while the former can keep growing inside the cosmic horizon. This means that, for higher 
£, the vector mode dominates the temperature anisotropics of the CMB over the scalar and tensor modes from the 
PMF and the contribution from primary anisotropics. The integrated amplitude of the gravitational waves from the 
PMF can be negligible after horizon crossing. This is because the homogeneous solution begins to oscillate inside the 
horizon and decay rapidly [H, [13, [H, [11] . Consequently, gravity waves only affect the anisotropy spectrum on scales 
larger than the horizon at recombination. The tensor mode from the PMF therefore decreases at higher £. Thus, 
the vector mode of the CMB polarization from the PMF dominates for higher £ [Panel (Id) of Fig. 2] [O, [llj [i3| In 
the primary spectrum for higher £ the absolute value of the EE mode from the PMF is relatively-small [Panel (Ic) of 
Fig. 2] compared to the TT mode [Panel (la) of Fig. 2]. Hence, even though the EE mode of the primary spectrum 
damps less than the TT mode, the EE mode remains much smaller than the TT mode at higher £ in the flnal power 
spectrum. For the TE mode [Panel (lb) of Fig. 2], the contribution from the PMF vector and scalar modes can 
be comparable to that from the primordial curvature perturbations for higher £. Except for a small dip region near 
£ = 40, the scalar mode is always small compared to the primary spectrum. 

Regarding the BB mode, we note that the BB mode signal described in this paper is due to magnetic-ficld-induccd 
CMB fluctuations (with the peak around £ ~ 2000 as in [ll|, [11]). We do not include magnctic-fleld related BB 
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polarization coming from the Faraday rotation effect (with the peak around £ > 15000) discussed in [lO|, |39|. In our 
model, the BB mode from the PMF can dominate for £ > 200 if Bx > 2.0 nG. A potential problem in attempting 
to detect this signal on such angular scales, therefore, is the contamination from gravitational lensing which converts 
the dominant EE power into the BB mode [s^. However, since we already know quite accurately what the spectrum 
of the lensing signal must be, we can subtract its power directly. After removing the foreground effect, the BB mode 
from the PMF effect dominates for higher £ even for PMF parameters allowed by the CMB temperature constraint 
[1 . IITI . |44| . Note, that the there is a change of scale between panels Ic and lb. The BB mode and EE modes are of 
comparable magnitude. 

In Panel 2 of Fig. 2 we depict the CMB temperature and polarization anisotropics in the presence of a PMF taking 
into account the correlations. Since we obtain the temperature and polarization anisotropies of the CMB with the 
PMF from isocurvature initial conditions, the phase of the CMB perturbation with the PMF is different by 7r/2 from 
those without the PMF (on the adiabatic initial condition). The first, third, and odd numbered peaks of the scalar 
mode of the CMB perturbations rise for a positive correlation between the PMF and the energy density perturbations, 
while they are suppressed for a negative correlation. These are compared in Panel 2 of Fig ure 2 with the observed 
power spectra (WMAP [H,!!^, ACBAR [13, CBI DASI BOOMERANG [11, and VSA Panels 

(2a) and (2b) of Fig. 2 show clearly that the power spectral index of the PMF, ne, is more effectively constrained from 
CMB observations for lower £ than those for higher £. The models with higher ?7b give better fits to the observations 
than those with lower ne for the lower £ regions of the TT and TE modes. Furthermore, there is no discrepancy at 
higher £ between observations and theories of the CMB polarization for models with a PMF [Panels (2c) and (2d) 
of Fig. 2]. In our previous work [H, [13 there was a problem from the strong degeneracy between ub and Bx. This 
degeneracy, however, is broken by the different effects of the PMF on the CMB power spectrum for lower and higher 



The scalar-mode CMB tcmperature-anisotropy power-spectrum shape (^-scaling for different ne at large scales - 
low £) agrees with the results of [2^ , while the E-polarization power spectrum shape does not follow semi-analytical 
estimate in that paper. This difference is caused by the fact that we have included the effects of reionization [5l|, [s^l 



which were neglected in [26|. If the universe is re-ionized at z^e, CMB photons are scattered by electrons and the 
polarization is generated again. Since re-ionization results in a new scattering surface at relatively short distances 
from us at (and a relatively recent era), the viewing angle of the polarization from re-ionization becomes large. Thus, 
reionization causes some power to shii^t to lower multipoles [ssl . . 

We also point out that the observed power spectrum of temperature fluctuations in the CMB is likely to depend on 
frequency [23] • Such dependence is theoretically expected to originate from foreground effects such as the Sunyaev- 
Zel'dovich effect at higher multipoles. In contrast, the effects of a PMF are frequency-independent because the PMF 
affects the primary CMB as a background. Therefore, the correlation between the PMF and other foreground effects 
should be weak. Because of this, one should be able to eventually distinguish the PMF from foreground effects by 
using more than two observational data sets at different frequencies. 

In summary, we have found that we can constrain more precisely the power law index tib and amplitude B\ of the 
PMF from all modes of CMB temperature and polarization anisotropies. The strong degeneracy of these parameters 
[l6l . ITtI is broken by the different effects of the PMF on the CMB power spectrum for lower and higher £. The scalar 
mode from the PMF can be a main source for lower £, while the vector mode can dominate for higher £ in the CMB 
temperature anisotropies. Furthermore, these calculations suggest that it is possible to place a limit on the correlation 
parameters s'^^^ for large negative values of the spectral index. For example , s^"^\ s*^^^ < for ne = —2.9 is ruled out 
from the effects of the TT mode on both the lowest and highest multipoles as shown in panels 2a and 2b of Figure 2. 
Such results may constrain models for the origin of the PMF, along with other PMF parameters. 
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If the order of an integrand function f{x) oc x" is more than -1, we can integrate f{x) numerically and analytically, for 
example, J [1/ ^/x)dx oc ^/x 

They, however, have only showed the result for a power spectral index ne = and they have not analyzed nor explained 
the features of the approximation errors for jib < —1.5 in particular. Thus, we can not compare the our analysis of errors 
from the previous approximation method for na < —1.5 with the result of Brown and Crittenden [46|. 
We do not consider the magnetic- field-related BB polarization coming from the Faraday rotation effect p^. [3^. 
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FIG. 2; CMB temperature and polarization anisotropies from the PMF. Panels (la), (lb), (Ic), and (Id) show TT, TE, EE, and 
BB modes, respectively, for models with Bx ~ 4.0nG and ub ~ —2.5 or —2.9 as labeled. Panels (2a-d) show a comparison of 
the computed total power spectrum with the observed CMB spectrum for B\ = 4.0 nG and various values of ub, s''^\ ands^'^^ 
as labeled. Plots show various ranges for: a) TT(2 < i < 1000), b) TT(400 < £ < 3000), c) TE(2 < £ < 4000), and d) 
BB(2 < £ < 3000) modes Curves in all Panels are theoretical lines as indicated in the legend box on the figure. Lines in (lb) 
and (2c) are plotted in the absolute value. Downward arrows for the error bars of Panels (2c) and (2d) indicate that the data 
points are positive and the lower error negative. 



